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Abstract 

^q' Following the technique of [1], we show that strictly static and strictly 

stationary solutions of the Einstein-Maxwell equations are analytic in har- 
monic coordinates. This holds whether or not the Maxwell field inherits 
the symmetry. 
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ir~i ■ 1 Introduction 

(N 

In a classic paper from thirty seven years ago, [1], Muller zum Hagen showed 
that static solutions of the Einstein vacuum equations are analytic in suitable (in 
fact harmonic) coordinates, at least where the staticity Killing vector is time- 
like (which we shall call strictly static). He commented, without going into 
details, that the same should be true for the Einstein equations with sources, 
given restrictions on the kinds of source. In a second article, [2], he showed that 
stationary solutions of the Einstein vacuum equations are also analytic in suit- 
able coordinates, where the Killing vector is time- like (say strictly stationary; 
in this article we shall only consider strictly static or strictly stationary solu- 
tions). It has been part of the folklore of General Relativity, though doesn't 
seem to have been published explicitly, that the Einstein-Maxwell equations 
with a Maxwell field inheriting the symmetry have this property. (Recall that 
the Maxwell field inherits the symmetry if CkF = where K is the Killing 
vector generating the symmetry and F is the electromagnetic field tensor). In 
this note, we shall show that this is true, and that it is also true if the Maxwell 
field does not inherit the symmetry (i.e. when CkF ^ 0). 

One would like to know whether the analyticity extends up to or even across 
any horizons, as this would assist black-hole classification theorems. Known 
static or stationary black-hole solutions all are analytic across the horizon but, 
aside from [3] for static vacuum solutions with non-degenerate horizons, one 
does not have any direct proofs that this must be so. 
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The basic technique in this article is to follow [1] and use the following 
theorem of Morrey [4] : 

Theorem 1.1 (Morrey) Given a system of second-order equations 

$ A (x a ,u B ,u B a ,u B ab ) = 

in an unknown u A for A = 1,2, . . . , N , where & A is analytic in all arguments 
and the system is elliptic in a domain D, then any (7 2+M solution in D is in 
fact analytic. 

Recall that the system is elliptic if the symbol is non-degenerate, that is to say, 
define 

L A (x a ,p a ) = Y^V a p h ^B$ A (x a ,u B ,u B a ,u B ab ) 
and then we require 

detL A (x a ,p a )^0, 

for all x a in D and all nonzero p a . 

The method of proof will be just as in [1], that is we shall cast the Einstein 
equations into a form in which Theorem 1.1 can be applied. For strictly static, 
inheriting Einstein-Maxwell this is straightforward, but for the non-inheriting 
case there is a small complication. One of the Maxwell equations is first-order 
and so does not fit into the pattern considered by Theorem 1.1. Consequently, 
we raise the order and obtain a system to which Theorem 1.1 does apply, and 
then recover the first-order equation via the Bianchi identities. In the strictly 
stationary case, we first assume that a in equation (9) below is constant - an 
example shows that this extra assumption is necessary - then we raise the order 
and obtain a system to which Theorem 1.1 does apply. This time we don't 
seem able to recover the first-order Maxwell equation but still we can argue as 
follows: any solution of the system including the first-order Maxwell equation 
satisfies the system with the order raised and therefore is analytic. 

The plan of the article is as follows. Material on the static case is in section 
2 and on the stationary case is in section 3. In sub-section 2.1, we review 
the static Einstein equations with a general stress-tensor. In sub-section 2.2, 
we review some fairly standard theory of the static inheriting Einstein-Maxwell 
equations as it may be found for example in [5], and note that the equations are 
already in a form to which Theorem 1.1 can be applied. Thus these solutions 
are analytic. Then in sub-section 2.3, and drawing on [6], we recall the theory 
of non-inheriting Einstein-Maxwell equations and see how the argument above, 
first raising the order and then recovering the first-order equation, goes through. 
The corresponding programme is carried through for stationary solutions in 
section 3. 

2 The strictly static case 
2.1 The Einstein equations 

In [1], Miiller zum Hagen first shows that comoving coordinates can be intro- 
duced: 
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Theorem 2.1 (Miiller zum Hagen) If a static space-time M is C 5 , with the 
Killing vector K being C 4 , and the metric C 3 , then there is a C 4 atlas with 

g = F(x k )dt 2 - h ij {x k )dx i dx j 

with F and hij both C 3 . 

We shall write a,b, . . . for space-time indices and i,j, k, . . . for spatial indices. 

In the standard way (see e.g. [5]), for a convenient form of Einstein's equa- 
tions, we set F = e 2U , which is why we insist on strict staticity, and introduce 
the rescaled metric: 

hij = 6 hij . 

The (three-dimensional) Ricci tensor transforms as 

Rij = R tJ - DiDjU + U :i Uj - hij{AU + \DU\ 2 ), 

where Di is the Levi-Civita derivative for the metric hij and A = h lJ DiDj is 
its Laplacian. The Einstein field equations can be written in terms of h^ as the 
system 

Rij = DiDjU + U t U J + P t a P^T ab -^Tg ab ) 

AU+\DU\ 2 = e- 2U K a K\T ab -^Tg ab ), 

where Pf = gf is the projection into surfaces of constant t and K a is the Killing 
vector. Then, in terms of the rescaled metric hij the Einstein equations become 

Rij = 2U t U J + (P?P!;-h tJ e- 4U K a K h )(T ab -±Tg ab ) (1) 
AU = e- 4U K a K b (T ab -±Tg ab ), (2) 

where A is the Laplacian for h^. 

Still following [1], we introduce harmonic coordinates x l , i = 1,2, 3 for hij 
satisfying: 

h I (fl / hrnn% : m) ,n = ^ 

where h = det(hij), then (dropping primes) 

Rij = ~^h hij^ mn + Qij, (3) 

where Qij does not contain second derivatives of h^ . 

The vacuum Einstein equations can be obtained from (l)-(2) by dropping 
all terms in T ab . In harmonic coordinates, with the aid of (3), it is easy to 
check that they are in the form suitable for application of Theorem 1.1 and so, 
following [1], we conclude that strictly static solutions of the vacuum equations 
are analytic in these coordinates, which by their definition are only available 
away from any horizons (where F = 0). We seek next to carry through the 
corresponding calculation for the Einstein-Maxwell equations. 
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2.2 Inheriting Einstein-Maxwell 

In this section, we review some standard and fairly familiar material on the 
inheriting Einstein-Maxwell equations. This is in order to extend the result 
of the previous section to this case, but also to prepare the way for the non- 
inheriting case in the next section. 

Given a Maxwell field F ab which inherits the symmetry, we introduce the 
electric and magnetic fields in the usual way as 

E a = F ab K b , B a =*F ab K b . (4) 

(though some authors including [6] may put a factor e~ u on the right in these 
definitions). Inheritance of the symmetry is expressed as 

C-kEo, = 0, CxBa = 0, 

and, as is very familiar, implies that E a and B a are both closed. However, for 
a static Killing vector, one may without loss of generality set B a = (since 
staticity implies that the Poynting vector vanishes, so that E a and B a are 
proportional; the Maxwell equations then force the function of proportionality 
to be a constant and then a duality rotation sets B a to zero). 

In a simply-connected region, closedness of E a implies that there exists a 
potential <j) with E a = 4> A . The Maxwell tensor is 

Fab = 2e~ 2U Ey a K b -\, 

and, since K a is a static Killing vector we also have 

V a K b = 2U la K b] . 

The remaining Maxwell equation is now 

V a E a = 2E a U„, 



or 



V a (e- 2U E a ) = 0. 



In terms of d> this is 



V a (e- 2C/ ^a) = 0. (5) 
The energy-momentum tensor is 

ab — 2b ac b b + -^r cd t g ab 

= -2e~ 2U E a E b + 2e~ 4U \E\ 2 K a K b - g ab e- 2U \E\ 2 , 

where \E\ 2 = —g ab E a E b which is positive with our conventions. Substitute into 
(1) and (2) to find 



Rij = 2U, l U J -2e- 2U ^^ J (6) 



AU = < 2 '7,'.'o,o., (7) 
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and the Maxwell equation (5) becomes 

A<f> = 2h ij Ui(t>j. (8) 

The strictly static, inheriting Einstein-Maxwell equations have been reduced to 
the system (6)-(8). These equations are quite familiar and may be found for 
example on p285 of [5]. In harmonic coordinates, the equations are second-order 
elliptic and then by Theorem 1.1 their solutions are analytic. 



2.3 Non-inheriting Einstein-Maxwell 

In this section, we shall follow the derivations in [6] but with a slightly different 
definition of E a , B a and W a from that reference. If the static symmetry is not 
inherited, then, as shown in [6], without loss of generality, the Maxwell tensor 
satisfies 

C K F ab = -a-k F ab , (9) 
for real, non-zero constant a. With E a and B a defined as before, this gives 

C-REa = —aB a ; CxB a = aE a . 

We may therefore introduce W a by 

E a = W a sin(at) ; B a = - W a cos (at), 

to find 

C K W a = 0. 

The Maxwell equations become 

v b ( e - 2U w h ) = (10) 

e ab cd K d V a W b = -aW c , (11) 

and, in terms of W a , the energy- momentum tensor is 

T ah = -2e- 2U W a W b + 2e- w \W\ 2 K a K b - e- 2U \W\ 2 g ab . (12) 

From (1) and (2), with (10), (11) and (12), we obtain the system 

Rij = 2UiUj - 2e- 2U WiWj (13) 
AU = e~ 2V WWiWj (14) 
e t jk DjW k = ae- 2U W l . (15) 

When we rescale, we assume that Wi does not change; thus we must be careful 
to distinguish h^Wj from h^Wj. For the rest of this section, indices are raised 
and lowered only with h. 

The problem with the system just obtained is that equation (15) for Wi is 
first-order, so that Theorem 1.1 does not immediately apply. To deal with this, 
we decompose Wi into closed and co-closed parts: 

W t = <f> ti + fo, (16) 
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where D^ 1 = 0. Then from the divergence of (15), 

A<f> = 2h ij U >i {<j>j+iP j ), (17) 

which is a suitable, second-order elliptic equation for <p. The rest of (15) can 
be written 

e 2U e i jk b j ^ k = a{<j )ji + ^ i ). (18) 
Apply e 2U i k mi D m to (18) and simplify to obtain 

A^ fc = (5 fc 5 m -5 m 5 fc )V' m + 2^ m 5 fe V' m -2/ l mn ^ m 5„V'fc + a 2 (0,fc + V'fc), (19) 

where indices are raised and lowered with h. The first term on the right sim- 
plifies with the Ricci identity 

(b k b m - b m b k )^ m = 

and now we have a system consisting of (13), (14), (17) and (19) to which The- 
orem 1.1 can be applied. This system will have analytic solutions in harmonic 
coordinates, which is sufficient for our purposes: non-inheriting, strictly static 
Einstein-Maxwell solutions necessarily satisfy the system including (19) and so 
are analytic. 

It is worth noting that we can in fact recover the first-order equation (15) 
for Wj. To accomplish this, introduce the notation: 

so that (18) and (19) are respectively 

L^Lfyj = a^i + ^j). (20) 

Now introduce Xj by 

Xj = L j k W k -aW j , 
so that (15) is just Xj = 0- From (20) we have 

L i j Xj + a Xi = (21) 

so that either Xj = 0, which is (15), or Xi is a non-zero solution of (21), which 
is just (15) with the sign switched on a. Now return to the system (13)-(14) 
and calculate the contracted Bianchi identity. This turns out to be 

D'iRij - ^Rh tj ) = 2e- 4U i j ik W iX k, (22) 

where we have used 

DtW* - 2U ti W i = 0, 

which is permissible since, given (16) this is just (17) which is part of the system 
already solved. 



6 



Given the (analytic) solution of the system, the left-hand-side of (22) must 
vanish. Therefore so must the right-hand-side and either Xj = 0) i n which case 
(15) is satisfied, or Xi is proportional to Wj, say \i = /Wi f°r some function /. 
In this second case we obtain 

= L i J Xj + aXi 

= L/ifWj) I «f\l, 

= f(Xi + aWi) + afW i + W j L t if 

= f(f + 2a)W i + W j L i j f, 

using (21) and the definition of Xi- Now contract with W l and use the definition 
of L { 3 to find 

/(/ + 2a)W i W i = 0. 

If Wi isn't zero then / = —2a (since by assumption for this part neither x% nor 
/ is zero) but now 

Xl := L^Wk - aW 3 = -2aW t , 

so that 

L+Wk + aWj = 0, 

and we recover (15) but with the sign switched on a. Since a enters the second- 
order system we are solving (specifically (19)) only through its square, we should 
expect to have this ambiguity in the first-order equation (15). 

3 The stationary case 

We now follow the corresponding chain of arguments to see that, in the strictly 
stationary case too, first vacuum, then inheriting Einstein-Maxwell, and then 
non-inheriting Einstein Maxwell solutions, subject to an extra assumption, are 
analytic in harmonic coordinates. 

3.1 The Einstein equations 

We begin with the result from [2] on comoving coordinates in stationary space- 
times corresponding to theorem 2.1: 

Theorem 3.1 (Miiller zum Hagen) If a stationary space-time M is C 5 , with 
the Killing vector K being C 4 , and the metric C 3 , then there is a C 4 atlas with 

g = F(dt + Aidx 1 ) 2 - hijdx l dx j 

with F, Ai and hij all C 3 . 

As before, to simplify the Einstein equations we put F = e 2U and hij = e 2U hij, 
and introduce 

uji = e u e i jk d j A k . (23) 
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Our definition of to differs slightly from the definition in [5]: if the definition 
there is cof 5 then Wj = e~ 2U tof s . Note that 

DiJ -JUi = 0. (24) 

There is more choice in the coordinate t in the stationary case than in the 
static case. To reduce it, and following [2], we impose the condition that t be a 
harmonic coordinate, and then 

g ab V a V b t = e~ u Di{e u h l] Aj) = 0. (25) 

The Einstein field equations become 

Rij = DiDjU + UiUj + ^(uJitOj - hijUJ k u) k ) 
+P t a P b (T ab - l -Tg ah ) 
AU+lDUf + h^iViUj = e- w K a K h (T ah - l -Tg ah ) 
e^ k D 3 {e 2U u: k ) = 2e 3U P?K b (T ab - l -Tg ab ) 

These can be checked against the form of the equations given in [5], taking 
account of the different conventions there. 

For stationary vacuum solutions, we set the stress-tensor terms to zero and 
write the equations in terms of the rescaled metric hij to find 

Rij = 2U ti Uj + ^oj i uj j (26) 

AU = ~h*juiUj (27) 
~e l lk (D j u Jk + 2U^ k ) = (28) 

(here uji and in the next set of equations Ai are both assumed not to change 
under rescaling) while (23), (24) and (25) become 

Wi = e^e^djA,, (29) 
h ij bi{e- 2U Uj) = (30) 
Di(h ij Aj) = 0. (31) 

In harmonic coordinates, the system consisting of (26) and (27) is in a suitable 
form to apply Theorem 1.1 but (28) isn't. We write it instead as an equation 
on Ai using (29) to find 

ti j (DiDjA k + 4Ui(DjA k - b k Aj) - (DiD k - D k Di)Aj) = 0, 

when, with the aid of the Ricci identity, it takes the right form to apply The- 
orem 1.1. Now, following [2], we conclude that strictly stationary vacuum 
solutions are analytic in harmonic coordinates. 
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3.2 Inheriting Einstein-Maxwell 

For inheriting Einstein-Maxwell, we introduce the electric and magnetic fields 
as in (4) and then both are closed, so that 

E a = F ab K b = V a a, B a = *F ab K b = V a /3, 

for two scalar potentials a and j3. In rescaled variables, the Einstein-Maxwell 
equations become 

Rij = 2UiUj + ^UiUj - 2 e - 2U (a i a j + ptfj) 

AU = ^(2e- 2U (a i a j +P, i p, j )-uj i u j ) 

~e i jk {D j u k + 2U, j u k ) = -Ae-^e^aj^ 
Aa = h ij (2U ii a J + Ui(3j) 
A/3 = WpUjPj-UiCtj), 

together with (29)-(31). 

As for vacuum, we eliminate cji in favour of Ai, then in harmonic coordinates 
the system is elliptic, Theorem 1.1 applies and solutions are analytic. 

3.3 Noninheriting Einstein-Maxwell 

In the stationary but non-static, non-inheriting case, we still have (9) but it is 
now an extra assumption that a be constant. (For this, see the discussion in [6]: 
it is possible to have a non-constant a with a null Maxwell field; there are pp- 
wave examples which are stationary, non-static and non-analytic in harmonic 
coordinates, for example the metric 

g = 2du(dv + A 2 ((du) - 2d(d(, 

with Maxwell field F = Ae l ^ u ^du A d( + c.c. for real constant A and non- 
constant but non- analytic f(u); the stationarity Killing vector is K = d/du 
and C K F = f * F so that a = -/'). 

We shall therefore assume that a is constant. We introduce E a and B a as 
before, though these are not now closed. We find 

E a + iB a = ( a e iat , 

where (, a , the counterpart of W a in section 2.3, is complex and CxCa = 0. We 
may decompose Ca into real and imaginary parts as 

Ca = Za + iria- 

There is a residual 'gauge freedom' under change of choice of time-coordinate: 
t' = t + fix'), C = e-^Ca, K = A a - / j0 . 
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The Maxwell equations split into divergence equations which can be written in 
the form 

AC - (Ui - iuJi + iaAi)C = 0, 

and curl equations: 

e i 3k ( D jCk - iaAjCk) = ae~ u d. 
In the rescaled variables, these are: 

A {h ij Cj ) - h ij {2Ui- icoi + iaAi) Q = (32) 

and 

~ei jk (Dj( k - iaAjCk) = ae~ 2U (33) 

where Q is assumed not to change under rescaling. It is straightforward to 
check that the divergence equations follow from the curl equations, and that 
both sets retain the residual gauge invariance. The field equations are the 
Einstein equations in the rescaled variables: 

Rij = 2U, l U J + ^oj l LO j -2e- 2U (Uj + ViVj) (34) 

AU = ^(2e- 2U (Uj + mVj)-^j) (35) 
l i jk {b j uj k + 2U, j u k ) = -4e- 2U ^ m (36) 

together with the Maxwell equations. We need to raise the order of the first- 
order Maxwell equations, (32) and (33), as before, so decompose Q into closed 
and co-closed pieces: 

d = v,i + Zi 

for D-co-closed Z^: DiZ 1 = 0. By imposing the divergence equations (32), we 
obtain the following elliptic equation for v: 

Au = h ij {2U :i + i{aAi - uJi)){vj + Zj). (37) 
We shall raise the order on the curl equations (33). Set 

A, = e 2U ~e l jk {D 3 Ck ~ iaAjCk), (38) 
then for the new second-order equation, we impose 

e 2C, eV fc (¥fc " iaAjkk) = a 2 Q. (39) 

which is a second-order equation for Zi, in which the highest-order derivative 
terms are 

e w {w m (bjbi - b t bj)z m - h^bjbiZrn). 

Here, as in (19), the Ricci identity allows the first part to be expressed without 
derivatives, while the second part is in the right form for Theorem 1.1. 

Note that (39) preserves the gauge- invariance and, if (33), holds then Aj = 
aQ and so (39) holds. Thus any solution of the system (32)-(36) also satisfies 
the system (34)-(37) together with (39). This latter system, with u>i eliminated 
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in favour of Ai, is second-order elliptic and so Theorem 1.1 can be applied. 
We conclude that strictly stationary, non-inheriting solutions of the Einstein- 
Maxwell equations with constant a are analytic in harmonic coordinates. 

It doesn't seem so straightforward in the stationary case as in the static to 
recover the first-order equation (33) (note that the ambiguity in the sign of a 
which was a feature of recovering the first-order equation in the static case, is 
not present in the stationary case since, for example, (37)- (39) all contain a to 
the first power). This difficulty seems to be connected to the fact that there 
are six real equations in (33) but one obtains only three real equations from the 
contracted Bianchi identity. 
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